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Abstract. We consider the numbers arising in the problem of normal ordering of 
expressions in boson creation and annihilation a operators ([a,a^] = 1). We treat a 
general form of a boson string (a^) r ™a s ™ ...(a' f ) r2 a S2 (a' ! ) ri a Sl which is shown to be associated 
with generalizations of Stirling and Bell numbers. The recurrence relations and closed-form 
expressions (Dobihski-type formulas) are obtained for these quantities by both algebraic and 
combinatorial methods. By extensive use of methods of combinatorial analysis we prove the 
equivalence of the aforementioned problem to the enumeration of special families of graphs. 
This link provides a combinatorial interpretation of the numbers arising in this normal 
ordering problem. 



1. Introduction 

In this paper we consider a pair of boson creation a) and annihilation a operators satisfying 
the commutation relation 

[0,0^ = 1. (1) 

These operators play a fundamental role in the formalism of second quantization in Quantum 
Mechanics and Quantum Field Theory (QFT) PQ , |2] , • Since the creation and annihilation 
operators do not commute serious problems with their ordering arise. A very convenient and 
well defined form of the operators depending on a and a) is the so called normally ordered 
form [4 . An operator is said to be in a normally ordered form if all creation operators stand 
to the left of the annihilation operators. The most important application field of the normal 
order is the QFT jHj. For a recent study of the interplay of the QFT, normal order and 
combinatorics see Ref.jS]. Procedure of normal ordering of the operator, i.e. moving all the 
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creation operators to the left with the use of relation Eq.fTj). is in general a nontrivial task. 
A first example is ordering of the power of the number operator [6 J [7 : 

n 

( a t a ) n = J2S{n,k){a^) k a k , (2) 

k=l 

where S(n, k) are the Stirling numbers of the second kind [Hj enumerating partitions of the 
set of n elements into k nonempty subsets, and satisfying the following recurrence relation 
S(n + 1, k) — S(n, k — 1) + kS(n, k) with initial values S(n, n) = S(n, 1) = 1. 
As the extension of this result we have considered operators in the form (a^) r a s (r, s -positive 
integers, r > s), for which a normally ordered form is given by 

ns 

[(a^) r a s ] n = (at)"^ S rAn, k){a ] ) k a\ (3) 

k=s 

where S r ^ s (n,k) are generalized Stirling numbers [Oj , [Ej , [H] , [12] > [E3 , [H| • This kind of 
formulas allow one to write the exponentials e^^ 0,3 in the normally ordered form and 
then easily calculate the coherent state expectation values which are of importance e.g. in 
Quantum Optics [I] . The clue of these calculations is the knowledge of the properties of the 
numbers S rtS (n, k). As they are of a combinatorial origin, the recurrence relations, Dobihski- 
type formulas, closed- form expresions and generating functions were extensively studied [TU] . 
In the following we further extend these results to normal ordering of a general boson string 
in the form (a)) rn a Sn . . . (a^) r2 a S2 (a^) ri a Sl , by establishing a link to special structures in 
enumerative combinatorics. This in turn gives us the rigorous demonstration of the properties 
of the generalized Stirling and Bell numbers arising in this problem. The construction 
of the graphs (the so called 'bugs') associated with these numbers provides a graphical 
interpretation of the normal ordering procedure. 



2. Generalized Bell and Stirling numbers 

In this section we define the generalization of ordinary Bell and Stirling numbers which arise 
in the solution of the normal ordering problem for a boson string. Given two sequences of 
positive integers r = (r l5 r 2 , . . . , r n ) and s = (s 1 ,s 2 , ■ ■ ■ , s n ) we let S TiS (k) be the positive 
integers appearing in the expansion 

Sl+S2 + --- + Sn 

(o t ) ^ ''o* n ...(o t ) ra a aa (o t ) ^1 a al = (o t )* , ^ S r , s (k)(a^) k a k , (4) 

k=si 

where d n = Y^=i( r i ~ s «)> "which we assume here to be non- negative. We observe that the 
whole theory can be carried through for d n negative, at the cost of minor adaptations, which 
however do not change the numbers involved. Note that the r.h.s. of Eq.(0J) is already 
normally ordered. 

We call S r ^ s {k) the generalized Stirling numbers of the second kind. The generalized Bell 
number is defined as the sum 

k=s± 
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In this notation the generalized Stirling numbers defined in Eq. (JHJ) correspond to a uniform 

n 

, r) and s = (s, s t . . . , s). 
We introduce the notation r l±l r n+1 = (r±, r 2 , . . . , r n , r n +i) and s l±l s n+ i = 
(si, s 2 , ■ ■ ■ , s n , s n+ i) and state the recurrence relation satisfied by generalized Stirling 
numbers S r)S (k) 

S r \Dr n+ i,sm n+ i(k) = / ] I t J (d n + k — j) Sn+ i-jS T ,s(k — j), (6) 

3=0 V 3 J 

where (l) p = I ■ (I — 1) • ... • (I — p + 1) is the falling factorial. 

One can give the derivation of Eq. (JHJ) by induction using the following consequence of Eq. (JTJ) 
(see the proof in (2J): 

^vy = E r ) (o P (« f )^ fe - p . (7) 

p=0 \ P J 

The full details of this approach can be consulted in [Tlj . 

Observe that the problem stated above can also be formulated in terms of the multiplication 
X and derivative D operators as they satisfy [-D,X] = 1. The representation of boson 
commutation relation with the D and X operators resembles the Bargmann representation 
[I], used in connection with coherent states. (Here we do not enter into that framework, with 
all the intricacies of the scalar product, hermiticity etc., as in our context only the algebraic 
properties matter.) Then Eq.(J3]) can be rewritten as: 

S!+S2 + ... + S n 

X rn D s " . . . X r2 D S2 X ri D Sl = X dn S r)S (k)X k D k . (8) 

k=si 

Acting with both sides of Eq.fjHJ) on the exponential function e x we get the identity 

X r "D Sn . . . X r2 D S2 X ri D Sl e x = x dn e x B TjS (x) (9) 

where 

Sl+S2 + ... + S n 

B r , s (x)= J2 S rMx k (10) 

k=si 

is the so called generalized Bell polynomial. Observe that the order of the so defined 
generalized Bell polynomial does not depend on r. Eq.fjHJ) gives the formula 

X d ^e x B rWrn+usWSn+1 (x) = X r ^D s ^e x x d "B r , s (x). (11) 

Using the well known commutation rule (equivalent to the Leibniz rule) D n e x f(x) = 
e x (D + I) n f(x) we get the recursive formula 

B r ur n+M+1 (x) = X s ^- d -{D + iy^X^B^x) (12) 

By taking coefficients of x k on both sides of Eq. (fT2"|) we obtain the recurrence relation for 
the generalized Stirling numbers of Eq. ((HJ) • 

Observe that the action of the l.h.s. of Eq.fjHJ) on e x may be calculated explicitly. To this 
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x 



end one first evaluates it on the monomial x n yielding YYj=i(dj-i + n ) 
turn easily gives the result of the action on the exponential function e x . 
With this observation, together with Eq.® we arrive at the extended Dobihski-type relation 
for generalized Bell polynomials 



B r 



m 



X 

ml 



(13) 



which by Cauchy's multiplication of series yields the expression for S T>s (k): 

1 k f k\ n 

= V E L (- 1 ) fe_m - II ' + ^-l)., (14) 
" m=0 \ / j=l 

An alternative, very similar demonstration of the above results can be carried through with 
the use of coherent states. These are defined for complex z, as \z) = e - ' 2 ' 2 / 2 Yl^Lo ^yl n )' 
where a^a\n) = n\n), a\z) = z\z) and (n\n') = 5 n y [4 . The |n)'s are called the number 
states. The coherent state matrix element of Eq.(JIJ) establishes a link to generalized Bell 
polynomials of Eq. (fTUj) : 



(z\(a 



^) r2 a S2 (a j ) ri a Sl 



(15) 



which after some algebra, provides an equivalent derivation of Eqs. (jl3j) and ()14|) . The 
first instance where the relation between the coherent state matrix elements and the Bell 
polynomials appears is Ref. [Zj, again for the generic case of Eq.(J2J), for which conventional 
Bell polynomials are obtained. 

We shall proceed now to give a combinatorial interpretation of the above results. The 
essence of subsequent paragraphs will be a graph-theoretical description of the problem. We 
define the structures (graphs) that are counted by the generalized Bell and Stirling numbers 
and then give a thorough combinatorial derivation of the recurrence relations, Dobihski- 
type formulas and other closed-form expressions. The Eqs. ((HJ), (fT3*j) and (|T4"|) will emerge 
from purely combinatorial considerations and this will permit the bijective identification of 
algebraic and combinatorial structures. 



3. Bugs, colonies, settlements and recurrence relations 

We introduce now a number of tools needed to describe the problem in the graph-theoretical 
language. 



Figure 1. A (3, 2)-bug 
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Definition 3.1 A bug of type (r, s) consists of a body and s legs. The body is formed by 
r linearly ordered empty cells. Each foot of the s legs is labelled with one number from an 
integer segment (m, m + s] := {m + 1, m + 2, . . . , m + s}, see FigO 

Definition 3.2 Consider a set of n bugs, the first one of type (n, si) and feet-labelled with 
labels in (0, si], the second of type (r 2 , s 2 ) with labels in (si, Si + s 2 ] and so on. A colony is 
one of the possible ways of organizing the bugs using the following procedure. The first bug 
has to stand over the ground. Once the (J — l)th bug is placed, the jth is placed by putting 
some (or none) of its Sj feet in the ground and each one of the rest in one of the empty cells of 
the bodies of the preceding bugs, see FigEl The pair of sequences (r, s), r = {r\, r 2 , . . . , r n ), 
s = (si, s 2 , ■ ■ ■ j s n)i carrying the information about the types of the bugs is called the type 
of the colony. The legs of the colony standing on the ground are called free. 
Assume now that there is a set of m empty cells in the ground. An m- settlement is a colony 
where each one of the feet corresponding to the free legs is placed in one of the ground 
cells. A surjective settlement is one where all the ground cells are occupied. The type of a 
settlement is defined to be the type of the subjacent colony. 



A 

1 

Figure 2. A colony of type (3,2,1,3;2,2,2,3) and 5 free legs. 
The main theorem of interest for us is: 

Theorem 3.1 The Stirling number S ryS (k), s\ < k < S\ + s 2 + . . . + s n , counts the number 
of colonies of type (r, s) having exactly k free legs. The Bell number B r s counts the number 
of colonies of type (r, s). 

Before proving it, we state the following 

Lemma 3.1 A colony of type (r, s) and with k free legs has exactly d n + k empty cells. 

Proof. The total number of cells of the colony is equal to Y^i=i r i- The number of occupied 
cells is equal to the total number of legs minus the number of free legs (J^ILi s « ~ ^ 
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Now we are ready to prove the Theorem 13.11 

Proof. Denote by C T>s (k) the number of colonies of type (r, s) with exactly k free legs. Since 
C(rv,si)(k) = S(r V>sl )(k) = 6(sx,k) it is enough to prove that the numbers C r>s (k) satisfy the 
same recursion as the generalized Stirling numbers of Eq. © • 

Crtt)r n +l,8hJ5n+l(^) ~ I „ + ) (d n + k — j) Sn+1 _j C VtS (k — j) (16) 



3=0 



J 



The l.h.s. is the number of colonies of type (r l±) r n+ i,s t±J s n +i) having exactly k free legs. 
We claim that in the right hand side the expression 



s n+l 
J 



(d n + k-j) Sn+1 „j C r , s (k-j) 



gives the number of such colonies where the (n + l)th bug has exactly j free legs. Obviously, 
this would prove the identity. We now prove our claim. In order to get a colony with k 
free legs the (n + l)th bug has to be placed in a colony of type (r, s) and k — j free legs. 
C r , s (k — j) is the number of such colonies. We choose the free legs of the (n + l)th bug 

in [ Sn+1 ) ways. Since by proposition (13.11) there are d n + k — j empty cells in the n-bugs 

V j J 

colony, (dn + k — j) Sn+1 -j gives the number of ways of distributing the rest of the feet of the 
(n + l)th bug into the empty cells. □ 




Figure 3. A 12-settlement of the colony in Fig.J2J) 

In the next section will follow a number of propositions clarifying the properties of structures 
in question. 



4. Counting settlements and Dobinski-type Relations 

We shall count now the number of m-settlements which will provide the link with Eqs.(|13J) 
and (pHj) viewed from the combinatorial perspective. 
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Theorem 4.1 Let p(m, r, s) be the number oj ' m- settlements of type (r, s). We have 

n 

p(m,r,s) = Y[(m + dj-i) Sj . (17) 

3=1 

Proof. There are {m) Sl ways of placing the feet of the first bug into the m ground cells. 
After placing the (j — l)th bug there are m + dj_i empty cells available (the previously 
placed bugs have provided YliZi r i empty cells and occupied YliZi s i cells). Then, there are 
(m + dj-i) s . ways of placing the Sj feet of the jth bug. □ 

Corollary 4.1 We have the polynomial identity 

n Sl+S2 + --. + Sn 

11(^ + ^-1)^= S *,*(k)(x)k- (18) 

3=1 k=si 

Proof. By the previous theorem, for an integer value of x the l.h.s. counts the number of 
x-settlements of type (r,s). S rtS (k)(x)k counts the number of ways of settling a colony of 
type (r, s) with k free legs in x ground cells. Then, the r.h.s. is another way of counting 
x-settlements. □ 

The exponential generating function of the surjective settlements is equal to the polynomials 

Sl+S 2 + ... + S n k S 1 +S2 + ... + S n 

B r , s {x)= £ S F)S {k)k\^= S rMx k . 
k=s\ k=s± 



Corollary 4.2 (Extended Dobihski-type relations) 
We have the identity 



oo n 



B r , s (x)e x = II( m + d i :)> -7- ( 19 ) 

m=s\ j=l 

Proof. Taking the coefficient of ^ of the left hand side we obtain 



k=0 \ J k=0 

By the previous corollary it is equal to the coefficient of 2-p in the right hand side. □ 
From Eq. (jl9|) we obtain 

oo n m 

B r , s (x) = e° II " + d >-i);h ( 2 °) 

m. 

m=s 1 j=l 

and 

oo ^ n 

S r>s = B T M = e- 1 Y,— ] + ^-0.j • (21) 



mi 

m=si j=l 
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Taking the coefficient of on both sides of equation (J2U)) we obtain the formula for the 
generalized Stirling numbers 

1 * ( h \ n 

= * (-ir m U( m+d ^u- (22) 

Evidently Eq.(j2DJ) is identical to Eq. ([IH|l and so are Eqs.(j22J) and (JHJ). This emphasizes again 
the already stated bijective correspondence between algebraic and combinatorial structures. 



5. Uniform colonies and settlements 



A colony or a settlement with all the bugs of the same type is called uniform. A uniform 
colony with n bugs of type (r, s) is called a colony of type (r, s) n . Following the notation of 
[TU] the corresponding Stirling and Bell numbers, enumerating uniform colonies of type 

n 

(r, s) n , are denoted respectively by S rtS (n,k) and B rtS (n). Clearly, for r = (r, r, . . . , r) 

n 

and s = (s,s t . . . ,s), S r:S (n,k) = S r!S (k) and B r ^ s {n) = B rs . The recursive formula, 
Dobihski-type relations and its consequences appearing here are natural extensions of those 
investigated in [HI] . 



Figure 4. A (3, l)-bug and corresponding small planar tree. 

The case s — 1 can be mapped into trees and forests. An (r, l)-bug can be identified with 
a planar tree, i.e. a tree where the leaves are all connected to the root and linearly ordered 
(see Fig0]). An increasing tree is one where the internal vertices are labelled with labels in a 
totally ordered set and the labels increase on any path from the root to any internal vertex. 
The uniform colonies with s = 1 corresponds to forests of increasing r-ary planar trees. The 
free legs are the roots of the trees (see FigEJ) • For r = 1, there is only one 1-ary increasing 
tree for each n. Then B\ t \{n) = B(n) is the ordinary Bell number. 

The exponential generating function of the r-ary planar increasing trees T r (x) satisfy 
the differential equation (see |15j . Chap. 5) y' = (y) r . From this we obtain T r (x) = 
— (r — Tjx, for r > 1. The generalized Bell number B r \(n) counts the number of 
r-forests with n internal vertices. By the exponential formula we obtain 

oo n 

V Br^n)^ = exp{ ^1 - (r - l)x - 1}. (23) 
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Figure 5. A uniform colony of type (3, l) 3 and the corresponding forest of increasing trees. 



We quote the explicit expression [TU] : 



T(n + 



r-l> 



(24) 



^r(i + ^)(fc-i)! 

In a subsequent publication we shall demonstrate that the summation formulas of the type 
Eq. (j24)l can be also obtained for many other strings describing the uniform case. 



6. Conclusions 



We have obtained analytic expressions and combinatorial interpretation of the integers 
generalizing conventional Bell and Stirling numbers, arising in the normal ordering of a boson 
string. All of their properties can be interpreted in terms of graph-theoretical language. The 
proof of the main result may also be obtained with the use of combinatorial theory of 
species [TH] , JT^J , [T7| . The results constitute an application of combinatorial analysis which 
produces the solution of quantum mechanical problem of normal ordering. For alternative 
interpretations of the numbers investigated in this work see Refs.JTK] and [THJ. It is an 
outstanding problem how to extend the key results of this work to the boson g-analogs. In 
this respect the Refs.|2Uj.|21j and will be of essential help. 
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